Coexistence of solutions in dynamical mean-field theory of the Mott transition 
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In this paper, I discuss the finite-temperature metal-insulator transition of the paramagnetic 
Hubbard model within dynamical mean-field theory. I show that coexisting solutions, the hallmark 
of such a transition, can be obtained in a consistent way both from Quantum Monte Carlo (QMC) 
simulations and from the Exact Diagonalization method. I pay special attention to discretization 
errors within QMC. These errors explain why it is difficult to obtain the solutions by QMC close to 
the boundaries of the coexistence region. 



PACS numbers: 71.27.-fa 

The Mott-Hubbard metal-insulator transition is one of 
the fundamental problems in the field of correlated elec- 
tronic systems. During the last few years, dynamical 
mean- field theory (DMFT) (l) § has emerged as an ap- 
propriate paradigm for this transition. Within DMFT, 
the paramagnetic half-filled Hubbard model on the Bethe 
lattice was found to undergo a first-order phase tran- 
sition at finite temperature T — The hallmark 
of such a transition is the coexistence of two solutions 
for the same value of the electronic interaction U with 
UciiP) < U < Uc2{P)- As is usual, I will label "insu- 
lating" the solution which can be continuously followed 
at low temperature from large values of the interaction 
down to Uci and "metallic" the one which exists from 
[/ = up to . 

Recent work |Q has challenged the above scenario. Us- 
ing the Quantum Monte Carlo (QMC) method, rather 
than the Exact Diagonalization (ED) approach 0|, the 
authors of ^ were unable to detect coexistence in the 
DMFT equations. This is surprising, as the two algo- 
rithms were compared extensively to each other (cf, e.g. 
ref. 1^, section VI, Appendix C). The aim of the present 
paper is to discuss the conflicting numerical approaches 
(for related work, cf ||^). I show that coexisting metal- 
lic and insulating solutions can indeed be obtained using 
both methods in a transparent and completely consistent 
way. I pay special attention to the problem of discretiza- 
tion errors. It will become clear that the QMC method 
needs to use very fine discretizations in order to obtain 
the insulating solution close to Uc^. Inside the coexis- 
tence region, the problem disappears. All the numerical 
work in this paper relies on programs which are publicly 
available 

Given the large number of articles already published 
on the subject, I will not repeat the standard definitions 
for the half-filled Hubbard model on the Bethe lattice 
within DMFT. These can be found e.g. in ref. |||. The 
important parameter of the model is the bandwidth D 
(cf. H); I quote interactions and temperatures in units 



of the bandwidth D/\/2. 

The point of departure of the present paper is ref. 
where coexisting solutions of the mean-field equations at 
inverse temperature [3D/^/2 — 100 were found (within 
ED) for Uc^ < U < with UcAPD/V^ = 100) - 
3.3D/V2 and C/^, if3D/V2 = 100) = 3.8D/V2. 

For this paper, I choose as a reference point the value 
U = 3.55D/V2, I3D/V2 ^ 100. At this point, both 
metallic and insulating solutions can be found easily by 
ED [0. These solutions are given by bath Green's func- 
tions GoiT) and impurity model Green's functions G{t) 
(cf for definitions). 

As we are well within the coexistence region, both 
solutions correspond to the bottoms of deep basins of 
attraction. I have, e.g., perturbed both the metallic 
bath Green's function, t/™'^*(T), as well as the insulat- 
ing one, t/™*(T), and fallen back into the corresponding 
solutions after a few iterations. This situation changes as 
we approach the boundary of the coexistence region. For 
smaller U, for example, the insulating solution becomes 
less and less attractive, until the basin of attraction van- 
ishes at [/ci • This scenario simply corresponds to the 
familiar free-energy landscape at a first-order transition. 

Both the ED and the QMC discretize some compo- 
nent of the DMFT equations. In ED, the discretization 
parameter is the number Ug of sites of the quantum im- 
purity model, while in QMC the number L of slices ap- 
pears, with At = (3/L (cf. ||]). I have computed the 
insulating bath Green's function Q^"^ (r) at the reference 
point for — 5, 6, 7. Each of the solutions (for fixed Hs) 
is fully converged, in addition, the convergence with the 
number of sites is excellent. I stress in passing that 
the sites in the ED algorithm are chosen in an optimal 
way, very similar to what is done in Gaussian integration 
|S|. Exponential convergence in Ug has been reported 
121. Given that ^^"^(r,^^ = 7) differs by less than 0.005 
from GoiT, ris = 6), I am led to the proposition that both 
G^^^{T,ns = 7) and G™^'^{T,ns = 7) are essentially exact. 
This assertion would have to be refuted for at least one 
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of the solutions (in fact, for the insulating one) if we were 
to agree with the authors of ref. . 

The following reasoning will lead us to a very in- 
teresting result: consider a single iteration of the self- 
consistency loop 



^ms,met ^ /->new I 'i\ 



If, contrary to our assertion, Q™^(t) or Q™^^(t) were not 
self-consistent, we should be able to detect differences, 
say, between Q™^ and Q^'^^ - even by other methods 
than ED, for example by QMC simulation. 

In this sense, I have discretized both t/™* and t/J"^* 
and used them as inputs for single self-consistency loops 
(as in eq. (|l|)) of the QMC algorithm for different values 
of Ar. From the discretized Green's functions G^'^(t), 
we then compute Q^^^ ir^ by inverse Fourier transfor- 
mation. Results for the metallic solution are shown in 
Fig. 1. We see that the finite- At effects are quite small. 
Clearly, g^^'^[T, Ar) ^ Q'^''\t) as At ^ 0. The metal- 
lic solution has not been contested. 




FIG. 1. ef''^(r,Ar) of the metallic solution, as com- 
puted by QMC for ArD/^ = 1,0.5 and g™'=*(r,n, = 7) 
(from below). Parameters of the reference point are 
^5D|^f2 = 100, U = 3.55i3/\/2. 

The situation becomes much more illuminating as we 
consider the insulating solution Q'^^^(t). In this case, I 
have computed G"^^(t) for Ar = 1,0.5,0.25. The results 
for Q^^^ (t\ again obtained by inverse Fourier transfor- 
mation, are shown in Fig. 2: it is evident that very large 
finitc-Ar effects are present. Notice that all the curves 
Q^^^ {t^ At) expose a plateau C/o ^ constant for r away 
from or (3. Such a plateau is characteristic of an insu- 
lating solution (cf. Fig. 1), and its value decreases with 
the gap of the single-particle density of states. We thus 
arrive at the crucial observation that the finite-Ar effects 
bias the self-consistency condition of the QMC algorithm 
towards the metallic solution. Nevertheless, we find back 
the ED solution as At 0. 




FIG. 2. C/?^'^(r, Ar) of the insulating solution, as com- 
puted by QMC for ArD/^/2 = 1, 0.5, 0.25 and GtriT, ris = 7) 
(from below) at the reference point. 

In terms of the free-energy landscape mentioned ear- 
lier, this simply means that the picture at finite Ar is 
tilted. The QMC iteration is dragged away from the in- 
sulating solution CJq"'* into the metal. 

Besides this effect, there seems to be no difference be- 
tween the iteration loop of ED and the one of QMC. 
Within a deep basin of attraction of the insulating so- 
lution, a small drag due to finite Ar effects should only 
lead to a shift of the solution, and its stability should 
be preserved. This is exactly what I have observed. I 
have run full iteration loops starting from Q^''^{t) for 
AtD/\/2 = 1 and Ar£)/>/2 = 0.5. In the first case, 
the simulation moved away from the insulation solution: 
After about 10 iterations, the metallic solution was ap- 
proximately recovered. In contrast, for AtZ?/\/2 = 0.5, 
the insulating solution is very clearly stable. A 40-day 
simulation of this single problem on a work station has 
obtained a very well-converged self-consistent solution for 
Goir, Ar) within QMC. This solution makes no more ref- 
erence to ED, but it of course resembles the curves shown 
in Fig. 2; the plateau value is ^o(t£'/%/2 = 50) 0.30, 
comparable to the one-shot solution at the same value of 
Ar. 

In my opinion, the present discussion of the discretiza- 
tion errors and both the one-step iteration and the ex- 
plicit self-consistent Monte Carlo solution leave little 
room but to accept the coexistence at finite tempera- 
ture. It is evident that the Monte Carlo simulation at 
the insulating solution has important finite-Ar effects. 
These errors modify the qualitative aspects of the solu- 
tion only close to the phase boundaries where the basins 
of attraction of the insulating solution are shallow, and 
small. Therefore, the discretization errors will be more 
pronounced close to , and also at higher temperature. 
This is what happened in the simulation of ref. [Q, and 
what led the authors to wrong conclusions. 
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